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Fuel-cell-based  auxiliary  power  units  can  help  to  reduce  fuel  consumption  and  emissions  in  transporta¬ 
tion.  For  this  application,  the  combination  of  solid  oxide  fuel  cells  (SOFCs)  with  upstream  fuel  processing 
by  autothermal  reforming  (ATR)  is  seen  as  a  highly  favorable  configuration.  Notwithstanding  the  neces¬ 
sity  to  improve  each  single  component,  an  optimized  architecture  of  the  fuel  cell  system  as  a  whole  must 
be  achieved.  To  enable  model-based  analyses,  a  system-level  approach  is  proposed  in  which  the  fuel 
cell  system  is  modeled  as  a  multi-stage  thermo-chemical  process  using  the  “flowsheeting”  environment 
PRO/II™.  Therein,  the  SOFC  stack  and  the  ATR  are  characterized  entirely  by  corresponding  thermody¬ 
namic  processes  together  with  global  performance  parameters.  The  developed  model  is  then  used  to 
achieve  an  optimal  system  layout  by  comparing  different  system  architectures.  A  system  with  anode  and 
cathode  off-gas  recycling  was  identified  to  have  the  highest  electric  system  efficiency.  Taking  this  system 
as  a  basis,  the  potential  for  further  performance  enhancement  was  evaluated  by  varying  four  parameters 
characterizing  different  system  components.  Using  methods  from  the  design  and  analysis  of  experiments, 
the  effects  of  these  parameters  and  of  their  interactions  were  quantified,  leading  to  an  overall  optimized 
system  with  encouraging  performance  data. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  expansion  of  onboard  convenience  services  such  as  enter¬ 
tainment,  navigation  or  air  conditioning  and  the  substitution  of 
mechanical  or  hydraulic  components  with  electric  units  give  rise 
to  an  increasing  demand  for  supplementary  electric  power  in  road 
vehicles,  ships  and  aircraft.  Moreover,  airplanes  and  long-haul 
trucks  that  require  electricity  during  parking  periods  face  public 
and  legal  pressure  to  reduce  local  emissions.  Therefore,  auxiliary 
power  units  (APUs)  are  gaining  importance  in  the  transportation 
sector.  APUs  efficiently  supply  electric  power  in  the  idle  state  and 
they  avoid  performance  losses  in  the  main  engine  of  propulsion 
during  normal  operation.  Particularly  when  they  are  based  on  fuel 
cell  (FC)  technology,  auxiliary  power  units  increase  operation  effi¬ 
ciency,  generate  near-zero  emissions  and  are  operated  with  low 
noise  [1-4]. 

In  contrast  to  batteries,  FC-APUs  are  not  constrained  by  lim¬ 
ited  operation  ranges,  provided  that  the  fuel  cell  is  operated  with 
fuel  that  is  available  onboard.  This  implies  on-site  fuel  processing 
for  which  typically  one  of  the  three  main  routes  is  applied:  partial 
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oxidation,  steam  reforming  or  autothermal  reforming.  The  partial 
oxidation  process  (POX)  is  advantageous  since  no  water  is  needed 
for  the  reforming  process.  However,  partial  oxidation  of  diesel  or  jet 
fuel  often  leads  to  carbon  deposition.  On  the  one  hand,  Lawrence 
and  Boltze  [5]  successfully  conducted  start-up,  load-change  and 
shut-down  tests  with  a  sub-kW  diesel  fueled  SOFC  system  includ¬ 
ing  a  POX  reformer.  On  the  other  hand,  stable  operation  could  only 
be  maintained  for  a  few  hours.  Lindermeir  et  al.  [6]  tested  an  SOFC- 
APU  containing  a  two-stage  POX  reformer  exhibiting  260  W  net 
electric  power.  However,  using  an  alkane  blend  as  fuel,  no  more 
than  9  h  of  successful  operation  could  be  demonstrated. 

The  main  advantage  of  steam  reforming  is  its  high  efficiency 
due  to  a  higher  hydrogen  yield  compared  to  POX.  However,  during 
the  steam  reforming  of  desulfurized  diesel,  Mengel  et  al.  noticed 
a  degradation  of  the  catalyst  due  to  soot  depositions  on  the  reac¬ 
tor  components  and  heavy  deposits  on  the  catalyst’s  inner  surface 
even  within  run  times  as  short  as  4  h  [7].  Ming  et  al.  [8]  presented  a 
steam  reformer  for  ultra-low  sulfur  diesel  (ULSD)  operated  at  low 
temperatures  ranging  from  490  °C  to  600  °C.  Reportedly,  conver¬ 
sion  rates  that  partly  exceeded  100%  related  to  the  Q  products 
were  obtained  for  310  h  on  stream.  More  accurate  values  below 
1 00%,  achievable  by  additionally  determining  the  residues,  have  not 
been  reported.  O’Connell  et  al.  [9]  implemented  an  oxidative  steam 
reforming  process  to  increase  long-term  stability  and  to  reduce  car¬ 
bon  depositions.  This  approach  is  similar  to  autothermal  reforming, 
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Nomenclature 

c  element  of  design  matrix,  factor  level  indicator 

e  relative  effect 

E  effect  (variable  units) 

E®  theoretical  electric  potential  based  on  lower  heating 
value  (V) 

F  Faraday  constant  (9.6485  x  1 04  A  s  mol-1 ) 

Ag°  change  in  Gibbs  free  energy  due  to  reaction  at  stan¬ 
dard  pressure  (J  kg-1 ) 

A hR  reaction  enthalpy 

Kwgs  equilibrium  constant  of  water-gas  shift  reaction 
LHV  lower  heating  value  (J  kg-1 ) 

N  total  number  (of  simulation  runs) 

N  molar  flow  (mol  s-1 ) 

Nc  number  of  carbon  atoms  per  molecule 

Pe  electric  power  (W) 

R  universal  gas  constant  (8.314 J  mol-1 1<-1 ) 

T  temperature  (I<) 

U  voltage  (V) 

Uf  fuel  utilization  factor 

y  molar  fraction 

Y  response  variable  (variable  units) 

Greek  symbols 

first-law  efficiency  (based  on  LHV) 
r]\oss  lumped  cell  voltage  loss  (V) 

X  normalized  air-to-fuel  ratio  (lambda) 

v  stoichiometric  coefficient 

Subscripts 

b  (after)  burner 

c  cathode 

i  species  i;  factor  i 

in  inflow 

j  counting  variable  for  simulation  runs 

k  index  of  response  variable 

m  modified 

N  Nernst- 

s  stack 

sys  system 

out  outflow 


however  at  a  much  lower  air-to-fuel  ratio  (AFR).  Based  on  a  micro- 
channel  technology,  they  obtained  over  98%  conversion  for  diesel 
during  38  h  of  operation. 

Autothermal  reforming  (ATR)  can  be  considered  as  a  compro¬ 
mise  between  partial  oxidation  and  steam  reforming.  Kang  et  al. 
[10,11]  experimentally  investigated  solid  oxide  fuel  cell  systems 
with  diesel  autothermal  reforming  reactors.  However,  fueling  the 
laboratory-sized  systems  with  a  binary  hydrocarbon  mixture  and 
with  sulfur-free  synthetic  diesel,  respectively,  resulted  in  rapid  and 
severe  degradation  due  to  incomplete  conversion.  Yoon  et  al.  [12] 
published  the  results  of  experiments  on  the  autothermal  reforming 
of  synthetic  diesel  fuel.  They  operated  their  reformer  for  approx¬ 
imately  lOOOh  with  an  LHV  efficiency  of  60%  and  observed  no 
significant  degradation  in  the  performance  of  the  reactor  over  time. 
Roychoudhury  et  al.  [13]  tested  their  microlith  ATR  reformer  with 
jet  fuel  (JP-8)  containing  low  sulfur  (~15ppm  S)  for  llOOh  on 
stream.  An  analysis  of  the  reformate  showed  that  propane  and 
propylene  products  remained  below  10  ppm,  resulting  in  a  con¬ 
version  rate  of  >  99.98%.  Long-term  experiments  on  autothermal 
reforming  with  commercial  premium  diesel  and  low-sulfur  jet  fuel 
(Jet  A-l,  <10ppmwS)  were  performed  by  Pasel  et  al.  [14]  and  Pors 


et  al.  [  1 5  ].  They  showed  a  conversion  of  99.7%  after  1 000  h  of  oper¬ 
ation  with  diesel  and  99%  after  2000 h  with  Jet  A-l.  In  order  to 
ensure  high  measurement  accuracy,  fuel  conversion  was  deter¬ 
mined  based  on  carbon  residues  (C2H4,  C2H6,  C3H7,  C3H8,  C4+Hx)  in 
the  dry  gas  phase  and  on  the  organic  carbon  content  of  the  conden¬ 
sate  withdrawn  from  the  cooled  product  gas.  At  1 1 00  h  of  operation 
with  Jet  A-l ,  the  mass  of  total  organic  carbon  (OC)  in  the  condensate 
was  analyzed  and  found  to  be  70  mg  OC  per  liter,  corresponding  to 
a  conversion  of  99.986%. 

Clearly,  the  best  reforming  performances  to  date  have  been 
achieved  using  autothermal  reforming.  In  comparison  to  POX 
and  SR,  longer  operation  times  and  higher  conversion  rates  were 
achieved  for  middle  distillates.  Therefore,  in  this  work,  autother¬ 
mal  reforming  was  selected  as  the  preferred  route  of  reforming. 
Furthermore,  different  fuel  cell  technologies  are  available,  ranging 
from  low-temperature  polymer  electrolyte  fuel  cells  to  high- 
temperature  solid  oxide  fuel  cells.  In  the  present  study,  we 
investigate  APU  systems  on  the  basis  of  solid  oxide  fuel  cells 
(SOFCs).  Despite  long  start-up  times,  SOFCs  have  several  advan¬ 
tages:  there  is  no  need  for  external  cooling  cycles,  complex  water 
management  is  avoided,  and  SOFCs  impose  the  least  stringent 
requirements  on  the  upstream  fuel  processing  compared  to  other 
fuel  cell  technologies  [16,17]. 

The  use  of  numerical  simulation  in  product  development  is 
becoming  increasingly  important,  since  prototyping  is  expensive, 
time-consuming  and  incapable  of  exploring  a  greater  number  of 
design  points  and  operating  conditions.  Consequently,  numerical 
modeling  of  fuel-cell-based  APU  systems  for  transportation  is  an 
active  field  of  research.  Petruzzi  et  al.  [18]  modeled  and  simu¬ 
lated  an  SOFC  system  for  APU  applications  but  excluded  the  fuel 
processing  steps.  Dobbs  et  al.  [19]  analyzed  an  SOFC-based  APU 
with  autothermal  reformer,  using  the  commercial  software  GCTool 
from  Argonne  National  Laboratory.  Lawrence  and  Boltze  [5]  numer¬ 
ically  examined  the  preliminary  design  of  an  SOFC-based  APU 
with  Matlab/Simulink®.  Matlab/Simulink®  was  likewise  used  by 
Sorrentino  and  Pianese  [20]  for  the  control-oriented  design  and 
analysis  of  a  hybrid  SOFC  APU.  Baek  et  al.  [21  ]  examined  SOFC  sys¬ 
tems  including  an  autothermal  reformer  based  on  a  Matlab®  model 
of  the  SOFC,  which  was  combined  with  an  Aspen  HYSYS®  represen¬ 
tation  of  the  system.  Baratto  et  al.  [2]  built  an  SOFC-APU  model  in 
AspenPlus®  as  a  basis  for  cost  estimation.  Most  recently,  Santarelli 
et  al.  [4]  conducted  flight  mission  simulations  of  SOFC  APUs  based 
on  a  zero-dimensional  system  model,  assuming  that  natural  gas 
was  available  as  fuel. 

Many  additional  numerical  studies  can  also  be  found  for  SOFC 
systems  based  on  natural  gas  or  biogas.  Schlitzberger  et  al.  [22]  and 
Farhad  et  al.  [23]  performed  simulations  of  SOFC  systems  based  on 
their  own  codes.  Bove  et  al.  [24]  developed  a  highly  parameterized 
black-box  SOFC  model,  which  may  be  included  in  standard  process 
simulation  software  such  as  AspenPlus®.  Palsson  [25]  integrated  a 
two-dimensional  finite-volume  model  of  the  stack  into  AspenPlus® 
to  enable  simulation  of  an  integrated  SOFC  and  gas  turbine  system. 
Zhang  et  al.  [26]  proposed  an  SOFC  model  based  solely  on  existing 
AspenPlus®  units,  which  nevertheless  requires  a  separate  calcula¬ 
tion  of  the  stack  voltage  [27].  Similar  models  were  developed  by 
Panopoulos  et  al.  [28]  and  Doherty  et  al.  [29]  in  order  to  simulate 
biomass-fueled  SOFC  systems.  A  study  including  off-design  anal¬ 
ysis  was  presented  by  Riensche  et  al.  [30].  The  authors  aimed  to 
optimize  a  200  kW  combined  heat  and  power  (CHP)  plant  using 
the  process  simulation  software  PRO/II™,  in  which  a  user-defined 
stack  model  was  included. 

Despite  the  variety  of  available  models,  all  of  these  approaches 
attempt  to  reproduce  the  characteristic  current-voltage  curves  of 
fuel  cells.  This  is  accomplished  by  postulating  reaction  kinetics  or 
equivalent  electric  circuits,  which  inevitably  bear  highly  param¬ 
eterized  correlations  that  depend  on  geometric  data  as  well  as 
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on  operation  characteristics.  The  latter  have  to  be  determined 
experimentally  for  each  particular  fuel  cell  stack.  However,  during 
preliminary  system  development,  final  stack  configurations  have 
yet  to  be  determined  and  there  is  often  no  measured  data  available. 
Moreover,  one  may  attempt  to  conduct  a  general  system  analy¬ 
sis  independent  of  a  particular  stack  but  focused  on  the  system 
components  and  their  interdependencies.  For  these  purposes,  we 
developed  a  new  SOFC  system  model,  which  is  presented  in  the 
following.  It  comprises  all  major  components  of  integrated  fuel  cell 
systems  and  allows  for  simple  and  generic  specification  of  the  SOFC 
stacks. 

2.  Calculation 

2.1.  Scope  and  modeling  principle 

The  auxiliary  power  units  under  consideration  consist  of  an 
SOFC  stack,  a  reformer,  an  anode  off-gas  burner  (afterburner),  a 
number  of  heat  exchangers  for  heat  integration  and  additional 
balance-of-plant  components.  These  various  system  components 
have  to  be  incorporated  into  a  single  simulation  framework  that 
must  be  simple  enough  to  limit  computational  effort  without  com¬ 
promising  comprehensiveness,  especially  with  regard  to  system 
interactions.  The  model  was  built  into  commercial  process  simula¬ 
tion  software,  so  that  existing  material  databases  and  pre-defined 
thermodynamic  property  calculation  procedures  could  be  used. 
PRO/II™  from  Invensys  Simsci-Esscor  [31]  was  used  for  this  pur¬ 
pose.  As  far  as  possible,  the  APU  model  was  implemented  based 
on  existing  unit  operation  modules  contained  in  the  software. 
This  enhances  design  flexibility  and  facilitates  comparative  studies 
including  variations  of  the  system’s  process  layout. 

2.2.  Solid  oxide  fuel  cell 

In  this  work,  the  operation  of  fuel  cells  with  reformate  is 
assumed,  so  that  hydrogen  and  carbon  monoxide  constitute  the 
effective  fuel  gases.  As  shown  by  Achenbach  and  Riensche  [32]  and 
Andreassi  et  al.  [33],  CO  is  directly  oxidized  at  SOFC  anodes.  Conse¬ 
quently,  two  net  electrochemical  reactions  have  to  be  considered: 

2H2  +  02^2H20  (1) 

2C0  +  02^2C02  (2) 

Furthermore,  due  to  the  catalytic  activity  of  the  Ni-containing 
anode,  the  following  two  equilibrium  reactions  proceed  simulta¬ 
neously  [34]. 

Steam  reforming  of  methane: 


Fig.  1.  Schematic  flow  sheet  of  the  SOFC  sub-model  including  Gibbs  minimization 
reactors  (Rl,  R2,  R3),  splitters  (SI,  S2,  S3),  controlled  gas  streams  (u,  e,  o)  and  heat 
fluxes  (dashed  lines). 

A4.  Effective  stack  voltage  and  overall  stack  current  adjust  inde¬ 
pendent  of  each  other. 

Assumption  A2  is  generally  in  agreement  with  the  literature 
(cf.  [27,32]).  Although  assumption  A3  is  in  contrast  to  the  fact  that 
methane  steam  reforming  is  kinetically  controlled  and  subject  to 
rate-limiting  diffusion  effects  [32,35],  it  nevertheless  provides  a 
satisfactory  approximation  [36,37]  and  automatically  considers  the 
possible  reverse  reaction,  i.e.  methanation.  As  will  be  shown  later, 
assumption  A4  implies  the  specification  of  an  effective  stack  voltage 
and  the  definition  of  a  fixed  fuel  utilization  factor. 

For  modeling  purposes,  the  SOFC  stack  is  conceived  as  a  multi¬ 
stage  process,  the  flow  sheet  of  which  is  outlined  in  Fig.  1.  First, 
reformate  and  air  are  brought  to  the  specified  operation  temper¬ 
ature  of  the  stack  according  to  Al.  Each  of  the  three  consecutive 
reaction  stages  is  represented  by  an  isothermal  Gibbs  minimization 
reactor ,  i.e.  the  equilibrium  composition  of  the  reaction  products  is 
calculated  based  on  the  principle  of  minimizing  the  Gibbs  energy 
at  constant  temperature  [38,39].  In  the  first  reactor  (denoted  ‘Rl  ’  in 
Fig.  1 ),  the  water-gas  shift  reaction  and  methane  steam  reforming 
simultaneously  proceed  to  equilibrium.  Subsequently,  incomplete 
fuel  utilization  is  accounted  for  by  dividing  the  fuel  gas  into  two 
branches  in  splitter  ‘SI’.  The  splitting  factor  is  equal  to  the  fuel 
utilization  factor  undefined  as  follows: 

_  NH2,out  +  NcO,out 

1  Nh2+Nco 


CH4  +  H20  -*  CO  +  3H2 

(3) 

Water-gas  shift  reaction: 

CO  +  H20  ->  co2  +  h2 

(4) 

In  the  present  model,  the  species  mass  balances  and  the  energy 
balance  of  the  complete  fuel  cell  stack  are  solved  as  a  function 
of  reaction  rates  and  reaction  enthalpies  based  on  the  reactions 
(l)-(4).  In  other  words,  the  SOFC  is  represented  by  a  single  con¬ 
trol  volume.  However,  the  progression  of  the  reactions  strongly 
depends  on  the  internal  geometric  structures  as  well  as  on  the 
material  properties.  Thus,  in  order  to  sustain  a  zero-dimensional 
model,  four  additional  assumptions  and  simplifications  are  intro¬ 
duced: 

Al.  Ideal  isothermal  operation  of  the  entire  fuel  cell  stack; 

A2.  Water-gas  shift  equilibrium  is  prevalent  everywhere  in  the 

stack; 

A3.  Methane  steam  reforming  proceeds  to  equilibrium; 


Once  specified  as  a  fixed  input  parameter  to  the  model,  Uf  further 
determines  the  amount  of  oxygen  that  permeates  through  the  solid 
electrolyte: 

^o2,o  =  2U/(^h2  +^co)  (6) 

As  can  be  seen  in  Fig.  1 ,  the  cross-electrolyte  flow  is  represented  by 
a  stream  of  pure  molecular  oxygen  (‘o’),  which  is  separated  from 
the  incoming  air  stream  in  stream  calculator  ‘S3’. 

In  the  fuel  cells,  the  oxidation  reactions  release  heat  and  electric 
power.  The  ratio  of  electric  power  to  the  total  energy  release  can  be 
described  by  a  first-law  efficiency,  which  depends  on  voltage  and 
operation  temperature  (Ts): 


<,0Ts)  = 


EUTs- 


i  =  H2 ,  CO 


(7) 


The  denominator  in  Eq.  (7)  contains  the  theoretical  electric  poten¬ 
tials  of  electrochemical  oxidation  of  hydrogen  or  carbon  monoxide, 
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respectively.  These  are  calculated  based  on  tabulated  data  of  terms  in  Eqs.  (lib)  and  (11c),  which  are  correlated  by  the  water-gas 
temperature-dependent  heat  of  formation  [40] :  shift  equilibrium: 


£g,,(Ts)=  AyTs),  1  =  H2,  CO  (8) 

In  the  present  model,  Eq.  (7)  is  used  to  define  a  stream  that  only 
contains  H2  and  CO.  This  stream,  denoted  ‘e’  in  Fig.  1,  is  removed 
from  the  utilized  fuel  gas  stream  ‘u’  and  is  subsequently  supplied  to 
reactor  ‘R2’  where  complete  oxidation  takes  place.  Consequently, 
the  resulting  reaction  enthalpy,  i.e.  the  heat  duty  of  ‘R2\  is  equal  to 
the  electric  power  delivered  by  the  stack: 

Pe  =  U-2¥uf(NH2+NC0) 

=  -UF[(^H,H2^H2)Ah^H2  +(*7h,c  o^co)Ahfl?co]  (9) 

The  remaining  fuel  gas  is  fed  to  reactor  ‘R3\  the  calculated  enthalpy 
difference  of  which  thus  corresponds  to  the  heat  released  inside  the 
stack.  A  reaction-independent  minimization  of  the  Gibbs  energy 
is  carried  out  [39],  so  that  the  product  gas  leaves  ‘R3’  in  overall 
chemical  equilibrium,  including  the  water-gas  shift  reaction. 

As  initially  stated,  all  reaction  stages  are  modeled  as  isother¬ 
mal  units.  On  the  other  hand,  the  stack  as  a  whole  is  assumed  to 
be  ideally  insulated.  Applying  the  energy  conservation  principle  to 
the  stack  consequently  imposes  a  constraint  on  the  only  remaining 
independent  variable,  which  is  the  flow  rate  of  air  to  the  cathode. 
The  air  flow  rate,  which  is  required  to  sustain  the  specified  outlet 
temperature,  is  expressed  in  terms  of  lambda,  i.e.  the  normalized 
air-to-fuel  ratio: 


Ac 


yo2Nair 

—  >  i 

Nq2,o 


(10) 


So  far,  it  was  assumed  that  a  constant  voltage  is  sufficient 
to  describe  the  fuel  cell  characteristics.  However,  if  the  fuel  gas 
compositions  change  due  to  alterations  of  the  process  layout,  this 
assumption  is  oversimplifying  and  will  thus  distort  comparative 
analyses.  Therefore,  a  voltage  correction  is  introduced,  which  is 
based  on  the  composition  dependency  of  the  Nernst  voltage.  Gener¬ 
ally,  the  actual  voltage  U  is  related  to  the  Nernst  voltage  UN  through 


U  =  UN-  1/1  ,oss  =  ^  In  (n^'j  -  ^loss  (Ha) 

The  definition  of  the  Nernst  voltage  contained  in  Eq.  (11a)  is 
based  on  the  assumption  of  ideal  gases,  isothermal  and  isobaric 
conditions.  The  change  of  the  specific  Gibbs  energy  is  a  reaction- 
dependent  constant,  while  the  logarithmic  term  accounts  for  the 
effect  of  the  actual  mixture  composition,  expressed  in  terms  of  the 
molar  fractions  y,  of  the  reacting  gases  and  the  stoichiometric  coef¬ 
ficients  Vi  of  the  corresponding  reaction.  The  lumped  voltage  drop 
rj\oss,  which  takes  into  account  the  ohmic  losses  and  the  cathodic 
and  anodic  overpotentials,  does  not  have  to  be  determined  explic¬ 
itly  within  the  present  modeling  framework.  Rather  it  is  important 
to  note  that  i)\oss  is  assumed  to  be  unaffected  by  composition 
changes.  With  regard  to  the  electrochemical  reactions  relevant  in 
the  SOFC,  i.e.  oxidation  of  H2  and  CO,  Eq.  (11a)  takes  the  following 
forms: 


u  = 

2  F 

*7loss,H2 

RT.  I 
2F  n  I 

(  yw2o  \ 

^yH2yg25 J 

(lib) 

u  = 

-Agc°o 

-  hlossXO  - 

2F  n  I 

(  yco2  \ 

(lie) 

2  F 
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Both  reactions  generate  the  same  voltage,  which  is  due  to  the 
equipotential  character  of  the  electrode  surfaces.  Thermodynami¬ 
cally,  the  voltage  equality  is  assured  by  the  composition-dependent 


f^WGS  = 


yco2yH2 

ycoyH2o 


(12) 


Eq.  (12)  contains  the  equilibrium  constant  JCW gs,  which  only 
depends  on  temperature  and  with  which  Eq.  ( 1 1  c)  can  be  rewritten: 


tt  _^co  RT  RT  (  yu2o 

U  =  — 2 p - dioss,co  ~  2 p  ln(^WGs)  ~  2 p  In 


yn2yo25 


(13) 


Comparing  Eqs.  (lib)  and  (13)  clearly  shows  that  if  the  composi¬ 
tions  of  the  anode  or  cathode  feed  gas  are  altered  as  a  result  of  a 
modified  process  layout,  the  voltages  resulting  from  both  reactions 
are  offset  uniformly.  The  corrected  voltage  Um  can  thus  be  deter¬ 
mined  as  a  function  of  the  specified  reference  voltage  U  and  the 
composition  alterations  of  H2,  02  and  H20: 


rr  rr  RT  1 

Um  =  U  -  —  In 
2  F 


(  ym, h2o  y^y°o2\ 

yym,H2ym502  yH2°  J 


(14) 


The  voltage  correction  accounts  for  deviations  from  a  reference 
state  that  is  associated  with  the  molar  fractions  yf  of  the  feed  gases, 
and  for  which  the  reference  voltage  U  is  defined.  If  the  feed  gas 
composition  is  changed  fromy*  toym  i-,  the  modified  voltage  Um  has 
to  be  considered  instead  of  U.  Strictly  speaking,  the  molar  fractions 
contained  in  Eq.  (14)  refer  to  local  conditions.  For  the  present  zero¬ 
dimensional  model,  however,  only  inlet  and  outlet  conditions  are 
known.  In  this  case,  taking  the  arithmetic  mean  values  provides  a 
satisfactory  approximation  [24]. 


2.3.  Reformer 

It  is  assumed  that  fuel  processing  is  accomplished  by  autother- 
mal  reforming  (ATR)  because,  as  stated  in  the  introduction,  it 
provides  the  best  long-term  stability  when  diesel  or  jet  fuel  is  used. 
Autothermal  reforming  comprises  both  exothermic  partial  oxida¬ 
tion 


CnHm  +  |o2^nCO+^H2 

(15) 

and  endothermic  steam  reforming 

C„Hm  +  n  H20  n  CO  +  (n  +  j)  H2 

(16) 

Furthermore,  side  reactions  have  to  be  taken  into  account,  most 
importantly  the  water-gas  shift  reaction  (reaction  (4))  and  metha- 
nation  [35,41]: 


CO  +  3H2  CH4  +  H20  (17) 

The  ATR  model  refers  to  an  apparatus  described  by  Pasel  et  al. 
[14],  albeit  without  an  integrated  steam  generator.  The  flow  sheet 
of  the  model  is  outlined  in  Fig.  2.  In  contrast  to  the  physical  appa¬ 
ratus,  air  and  liquid  fuel  are  mixed  first  before  steam  is  added. 
This  allows  the  internal  heat  transfer  from  the  reaction  zone  to  the 
two  cold  inlet  streams  (air  and  fuel)  to  be  modeled  using  a  single 
heat  exchanger.  The  heat  duty  of  this  unit  is  controlled  such  that 
the  temperature  of  the  gaseous  fuel-steam-air  mixture  is  equal  to 
measured  data. 

Since  exothermic  partial  oxidation  proceeds  much  faster  than 
endothermic  steam  reforming,  a  non-uniform  axial  temperature 
distribution  along  the  monolith  is  observed,  affecting  the  equilib¬ 
rium  reactions  (4)  and  (17)  [14,42].  Therefore,  the  ATR  sub-model 
comprises  two  reaction  stages  in  order  to  account  for  the  differ¬ 
ent  temperature  levels.  Each  reaction  zone  is  modeled  on  the  basis 
of  a  Gibbs  minimization  reactor.  In  the  first  reactor  unit,  denoted 
‘Zl’  (Fig.  2),  equilibrium  of  all  reactions  is  calculated  by  minimizing 
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Fig.  2.  Schematic  flow  sheet  of  the  ATR  sub-model  including  two  reaction  zones 
represented  by  Gibbs  minimization  reactors  (Zl,  Z2),  internal  heat  fluxes  (dashed 
lines)  and  control  of  heat  balance  (TC,  dotted  line). 

the  Gibbs  energy  and  assuming  adiabatic  conditions.  This  leads  to 
complete  conversion  of  the  liquid  fuel  due  to  partial  oxidation  and 
steam  reforming.  The  second  reactor  (‘Z2’)  is  an  isothermal  unit, 
whose  operating  temperature  corresponds  to  the  outlet  tempera¬ 
ture  of  the  autothermal  reforming  reactor.  The  temperature  drop  is 
accomplished  in  an  intermediate  heat  exchanger,  the  outlet  tem¬ 
perature  of  which  is  adjusted  by  means  of  a  controller.  The  latter 
is  specified  such  that  the  overall  heat  balance  of  the  ATR,  which  is 
assumed  to  be  ideally  insulated,  is  satisfied. 

2.4.  Auxiliary  system  components 

The  SOFC  stack  and  autothermal  reformer  are  the  two  main 
components  of  the  APU  system.  Nevertheless,  additional  system 
components  are  required  to  operate  the  fuel  cell  system.  The  fol¬ 
lowing  auxiliary  components  are  included: 

-  afterburner, 

-  heat  exchangers, 

-  condenser  and  separator, 

-  pumps,  blowers  and  compressors. 

An  afterburner  is  needed  to  ensure  low-to-zero-emission 
exhaust  since  the  fuel  gas  is  only  incompletely  utilized  by  the  fuel 
cells.  A  single  adiabatic  Gibbs  minimization  reactor  is  used  to  model 
the  burner.  The  calculated  product  temperature  corresponds  to  a 
theoretical  peak  temperature,  which  in  practice  should  be  kept 
below  950  °C  due  to  material  constraints.  This  is  accomplished 
by  controlling  the  amount  of  cold  excess  air.  The  corresponding 
lambda  value  is  defined  as  follows: 

Xb  =  - yof/air - ^ -  >1  (18) 

0.5(yH2  +ycO  +  4ycH4  )^anode-off 

Heat  exchangers  are  needed  mainly  for  the  purpose  of  heat  inte¬ 
gration.  Moreover,  heat  exchanger  units  have  been  introduced  in 
the  sub-models  of  Figs.  1  and  2  in  order  to  represent  internal  heat 
transfer.  To  all  of  these  units  an  overall  heat  transfer  coefficient 
of  100Wm-2I<-1  is  assigned.  With  the  use  of  micro-structured 
heat  exchangers,  even  higher  overall  heat  transfer  coefficients 
were  achieved  experimentally  [43].  At  150Wm_2K_1,  however, 
the  pressure  drop  already  exceeded  180  mbar.  Therefore,  the  more 
conservative  value  of  1 00  W  m-2 1<-1  was  assumed  in  order  to  com¬ 
ply  with  an  allowable  pressure  drop  of  50  mbar  (cf.  Table  1 ),  which 
was  reached  in  [43]  using  a  micro-structured  heat  exchanger.  All 
heat  exchangers  are  modeled  using  the  PRO/II™  unit  simple  heat 
exchanger.  In  order  to  be  able  to  estimate  the  required  minimum 
heat  exchanging  area,  ideal  countercurrent  operation  is  presumed. 

Autothermal  reforming  requires  a  constant  supply  of  water. 
Therefore,  the  APU  system  must  contain  a  closed  water  cycle  in 
order  to  be  self-sufficient,  which  is  realized  by  installing  a  com¬ 


Table  1 

Pressure  losses  assigned  to  individual  system  components. 


System  component 

Pressure 
loss  (mbar) 

System  component 

Pressure 
loss  (mbar) 

SOFC  stack 

Afterburner 

Anode  side 

20 

Combustion  chamber 

20 

Cathode  side 

30 

Heat  exchanger 

ATR 

Vaporization  of  water 

500 

Monolith 

70 

All  other  streams 

50 

Air  supply 

30 

Condenser 

Liquid  fuel  atomizer 

5000 

Moist  gas 

50 

Air  cooler 

3 

bined  condenser  and  water  separator.  This  water  recovery  unit 
is  modeled  on  the  basis  of  an  isothermal  flash ,  which  calculates 
the  vapor-liquid  equilibrium  at  constant  temperature  based  on  the 
non-random  two-liquid  (NRTL)  activity  coefficient  model  [44].  The 
solubility  of  supercritical  gases  is  considered  by  applying  Henry’s 
law.  The  temperature  of  the  flash,  which  is  needed  to  match  the 
mass  flow  of  condensed  water  with  the  water  requirement  of  the 
ATR,  is  adjusted  by  a  controller. 

Finally,  pumps,  blowers  and  compressors  are  needed  to  over¬ 
come  pressure  losses  in  the  system.  Blowers  and  compressors 
are  both  modeled  using  the  PRO/II™  unit  compressor,  to  which 
a  constant  isentropic  efficiency  is  assigned.  The  compressor  effi¬ 
ciencies  assumed  in  the  literature  typically  range  from  60%  to 
75%  [2,20,23,28].  In  this  work  an  isentropic  efficiency  of  60%  is 
used  for  the  base  case.  Since  this  is  a  rather  conservative  estimate, 
the  effect  of  increased  compressor  efficiency  on  the  overall  sys¬ 
tem  performance  will  furthermore  be  examined  in  section  4.2.  The 
efficiency  of  liquid  pumps  is  assumed  to  be  70%,  a  value  character¬ 
izing  well  dimensioned  equipment  [45].  The  back  pressures  that 
must  be  overcome  are  different  at  each  stream  inlet.  In  the  present 
study,  the  corresponding  values  are  determined  by  assigning  spe¬ 
cific  pressure  losses  to  all  flow  sheet  units  according  to  Table  1. 

3.  System  layouts 

The  first  part  of  our  study  comprises  a  comparative  study  of  dif¬ 
ferent  system  layouts.  A  base  system  is  built  as  a  reference  point. 
This  base  system  and  four  selected  alternative  designs  are  pre¬ 
sented  in  this  section.  They  will  serve  as  a  basis  for  the  process 
analysis  and  optimization  described  in  Section  4. 

3.1.  Base  system 

The  flow  sheet  of  the  reference  system  is  shown  in  Fig.  3.  The 
autothermal  reforming  reactor  is  fed  with  air  and  liquid  fuel  at 
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SOFC 


ambient  temperature,  while  water  enters  the  ATR  in  the  form 
of  superheated  steam,  which  is  generated  using  the  hot  burner 
exhaust  gas.  The  product  gas  of  the  ATR  (reformate)  feeds  the  anode 
of  the  SOFC  The  cathode  is  provided  with  moist  air,  which  is  pre¬ 
heated  by  hot  cathode  off-gas.  Ambient  air  is  assumed  to  contain 
20%  oxygen,  2%  water  vapor  (molar  basis)  and  nitrogen.  The  anode 
off-gas  is  directed  to  the  afterburner,  which  is  operated  with  addi¬ 
tional  ambient  air.  Fresh  air  is  used  in  the  afterburner  for  two 
reasons:  First,  the  burner  is  positioned  far  downstream,  so  that  a 
low-power  blower  can  be  applied  economically.  Second,  cold  air 
most  efficiently  limits  the  adiabatic  combustion  temperature  inside 
the  burner.  After  passing  through  the  heat  exchanger  to  generate 
superheated  steam,  the  exhaust  is  cooled  down  in  an  air-cooled 
condenser. 


In  contrast,  system  utilization  is  determined  by  a  control  volume 
that  encompasses  the  stack  as  well  as  the  recycle  loop. 

In  this  work,  three  configurations  including  recycles  are  exam¬ 
ined: 

-  Single  anode  recycle  (alternative  2). 

-  Single  cathode  recycle  (alternative  3). 

-  Combined  anode  and  cathode  recycle  (alternative  4). 

The  flow  sheet  of  alternative  4  is  shown  in  Fig.  5. 

4.  Results  and  discussion 


3.2.  System  variations 


4A.  Conceptual  study 


3.2 A.  Modified  operation  of  the  afterburner 

The  first  variation  comprises  the  operation  of  the  afterburner 
with  cathode  off-gas  instead  of  fresh  air.  Since  considerable  excess 
air  is  used  for  stack  cooling,  a  large  stream  of  cathode  off-gas  is 
available,  which  always  contains  enough  oxygen.  The  large  flow 
rate,  however,  results  in  a  significant  reduction  in  the  partial  pres¬ 
sure  of  the  water  vapor  in  the  burner  exhaust  gas.  Consequently,  the 
condensation  of  water,  which  is  necessary  to  close  the  water  cycle, 
becomes  challenging.  Thus,  the  condenser  is  installed  between  the 
SOFC  stack  and  the  afterburner  in  order  to  take  advantage  of  the 
higher  water  content  in  the  anode  off-gas.  The  resulting  system 
layout  of  alternative  1  is  depicted  in  Fig.  4. 

3.2.2.  Systems  including  recycles 

Three  more  system  variations  are  considered,  all  of  which 
include  recycle  streams.  Cathode  off-gas  recycling  can  generally 
be  used  to  replace  the  cathode-air  pre-heater  [30,46].  In  return,  a 
hot-temperature  blower  has  to  be  installed  in  order  to  overcome 
the  pressure  drop  inside  the  recirculation  loop.  The  recycle  ratio  is 
controlled  such  that  the  temperature  at  the  cathode  inlet  does  not 
fall  below  a  threshold  value,  which  is  determined  by  the  maximum 
allowable  difference  to  the  operation  temperature  of  the  stack: 
AT  =150  K. 

Recirculation  of  anode  off-gas  to  the  anode  inlet  is  applied  to 
reuse  a  portion  of  the  initially  unutilized  fuel  [21].  In  the  present 
study,  the  recirculation  ratio  of  anode  recycling  is  fixed  to  a  value 
of  0.5  -  a  tradeoff  between  increased  overall  fuel  utilization,  on 
the  one  hand,  and  decreased  stack  voltage  (Eq.  (14))  as  well  as  the 
increased  power  need  of  the  hot-gas  blower,  on  the  other  hand. 
In  the  presence  of  anode  recycling,  two  different  fuel  utilization 
factors  are  distinguished:  Stack  utilization  is  defined  by  the  inlet 
and  outlet  streams,  which  are  directly  attached  to  the  SOFC  stack. 


All  four  system  variations  and  the  base  system  were  evaluated 
separately  and  then  compared  to  each  other.  The  key  parameter  for 
evaluating  system  performance  was  the  net  electric  efficiency.  It  is 
defined  based  on  the  lower  heating  value  (LHV)  of  the  liquid  fuel: 

_  Pe, AC, net  _  ^DC/AC  '  ^e,DC  “  ^e,loss 

^sys  “  mfuel  •  LHVfuel  “  mfuei  •  LHVfuel  1  j 

The  parasitic  power  loss  Peloss  considers  all  electric  consumers  in 
the  system,  i.e.  pumps  and  compressors.  It  was  presumed  that  all 
system-intrinsic  consumers  run  on  alternating  current.  The  effi¬ 
ciency  of  the  DC/AC  converter  ^dc/ac  was  assumed  to  be  equal  to 
95%. 

For  the  evaluation  and  comparison  of  the  different  system  lay¬ 
outs,  equal  operating  parameters  as  well  as  operating  conditions 
were  used:  The  feed  streams  are  available  at  ambient  conditions, 
i.e.  at  1.01  bar  and  25  °C.  The  relative  humidity  of  the  surrounding 
air  was  assumed  to  be  60%.  All  systems  were  operated  with  a  fixed 
fuel  feed  of  1 .2  kg  n- dodecane  per  hour.  The  air  feed  to  the  ATR  is 
specified  by  the  ratio  of  molecular  oxygen  to  the  number  of  carbon 
atoms  contained  in  the  fuel: 


n(02)  _  Nq2 
"(C)  "Nfei’Nc 


(20) 


The  required  feed  of  water  is  defined  by  the  steam-to-carbon  ratio: 


n(H2Q)  Nh2o 
n(C)  JVfuel  ■  Nc 

The  values  were  set  to  n(O2)/n(C)  =  0.47  and  n(H20)/n(C)  =  1.9, 
respectively.  Under  these  conditions,  a  long-term  experiment  [14] 
demonstrated  the  stability  of  the  reforming  activity  with  99%  con¬ 
version  of  desulfurized  Jet  A-l  fuel  after  2000  h.  The  steam  further 
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Table  2 

Summary  of  simulation  results  of  the  conceptual  study. 


Base  case 

Alt.  1 

Alt.  2 

Alt.  3 

Alt.  4 

System 

Net  electric  efficiency  (%) 

29.2 

28.9 

31.4 

29.6 

32.0 

Net  electric  power  (AC)  (kW) 

4.30 

4.25 

4.62 

4.35 

4.71 

Parasitic  losses  (kW) 

0.54 

0.58 

0.65 

0.41 

0.49 

Fuel  utilization  (%) 

75.0 

75.0 

85.7 

75.0 

85.7 

SOFC 

Voltage  (mV) 

700 

700 

667 

691 

658 

Gross  elect,  power  (DC)  (kW) 

5.1 

5.1 

5.5 

5.0 

5.5 

Lambda 

8.7 

8.7 

9.4 

2.1 

2.2 

Air  standard  flow  rate  (m3  h-1 ) 

66.3 

66.5 

82.2 

15.8 

19.6 

Burner:  air  standard  flow  rate  (m3  h-1 ) 

5.9 

- 

2.0 

5.9 

2.0 

Overall  off-gas 

Standard  flow  rate  (m3  h-1 ) 

77.7 

72.0 

89.7 

27.2 

27.0 

Temperature  (°C) 

171 

201 

182 

417 

506 

Condensation  temperature  (°C) 

57 

59 

64 

57 

64 

Total  heat  exchanging  area  (m2) 

1.77 

1.88 

1.86 

0.86 

0.72 

has  to  be  supplied  to  the  ATR  at  a  specified  temperature  of  480  °C. 
The  reference  voltage  of  the  SOFC  stack  according  to  the  conditions 
in  the  base  system  was  defined  as  U = 0.7  V.  The  fuel  utilization  fac¬ 
tor  of  the  stack  was  always  equal  to  0.75,  and  the  uniform  operating 
temperature  was  assumed  to  be  750  °C. 

The  main  results  of  the  comparative  study  are  summarized  in 
Table  2.  Net  electric  efficiency,  net  electric  power  output  of  the 
system,  parasitic  losses  due  to  system  consumers  and  gross  elec¬ 
tric  power  of  the  SOFC  stack  correspond  to  Eq.  (19).  The  effective 
stack  voltage  was  calculated  using  Eq.  (14).  The  required  standard 
flow  rate  of  cathode  air  and  the  normalized  air-to-fuel  ratio  lambda 
are  related  to  each  other  through  Eq.  (10).  The  effective  off-gas 
comprises  the  exhaust  gas  from  the  afterburner  and  spent  cathode 
air.  The  indicated  condensation  temperature  is  a  threshold  value 
for  water  autarky.  Finally,  the  total  heat  exchanger  area  given  in 
Table  2  comprises  heat  exchangers  for  steam  generation,  conden¬ 
sation  and,  except  for  process  alternatives  3  and  4,  cathode  air 
pre-heating. 

Examining  alternative  1,  it  was  found  that  system  efficiency 
decreased  compared  to  the  base  case.  This  is  mainly  due  to  an 
increased  back  pressure  at  the  cathode  feed  line  as  a  consequence 
of  the  additional  pressure  drop  inside  the  afterburner.  Further¬ 
more,  the  first  alternative  design  requires  a  somewhat  larger  heat 
exchanging  area  due  to  smaller  temperature  differences  in  the 
steam  generator  and  air  pre-heater.  On  the  other  hand,  only  two 
instead  of  three  air  blowers  are  needed,  thus  reducing  costs  and 
space  requirements. 

The  net  electric  efficiency  of  alternative  2  increased  compared 
to  the  base  case.  Thus,  the  positive  effect  of  higher  system  fuel 
utilization  dominates  over  the  negative  effects  of  lower  stack  volt¬ 
age  and  increased  system-intrinsic  power  consumption.  Due  to  the 
increased  fuel  utilization,  more  excess  air  is  needed  to  hold  the 
stack  at  the  specified  operation  temperature.  On  the  other  hand, 
less  air  has  to  be  supplied  to  the  afterburner  since  the  anode  off¬ 
gas  contains  a  smaller  amount  of  combustible  compounds.  As  a 
consequence,  condensate  can  be  captured  with  less  effort  at  higher 
temperature. 

Cathode  off-gas  recirculation  as  realized  in  process  alternative 
3  leads  to  only  a  small  increase  in  the  total  electric  efficiency  com¬ 
pared  to  the  reference  case,  since  the  reduction  in  system-intrinsic 
power  consumption  is  largely  compensated  by  the  reduced  power 
output  of  the  stack.  The  latter  is  a  consequence  of  a  decreased  volt¬ 
age  due  to  oxygen  depletion  in  the  cathode  feed,  which  is  caused  by 
recycling.  The  flow  rate  of  fresh  air  to  the  cathode  was  reduced  by  a 
factor  of  four.  This  is  a  direct  consequence  of  the  high  recirculation 
ratio  (79%),  which  is  required  to  maintain  the  specified  opera¬ 
tion  temperature  of  the  stack.  Another  consequence  of  cathode 


recycling  is  that  the  off-gas  leaves  the  system  with  a  significantly 
reduced  flow  rate  but  at  elevated  temperature. 

Looking  at  system  alternative  4,  it  was  found  that  the  effects 
of  the  single  recycles,  which  were  examined  in  alternatives  2  and 
3,  approximately  add  up.  Compared  to  alternative  3  with  cathode 
off-gas  recycling,  the  lambda  value,  which  is  needed  for  isothermal 
operation  of  the  SOFC,  almost  remains  unchanged.  Nevertheless,  a 
larger  absolute  flow  rate  of  air  is  needed  since  the  fuel  utilization 
increases  due  to  the  anode  off-gas  recycling.  In  relation  to  the  base 
case,  the  net  electric  efficiency  increases  by  almost  10%,  while  the 
required  heat  exchanging  area  is  reduced  by  more  than  a  factor 
of  two.  The  latter  is  mainly  due  to  the  replacement  of  the  cath¬ 
ode  air  pre-heater  by  anode  off-gas  recycling,  which  can  neatly 
be  illustrated  in  a  temperature-enthalpy  diagram  using  compos¬ 
ite  curves.  Composite  curves  originate  from  pinch  analysis  with 
the  aim  of  identifying  the  maximum  possible  heat  recovery  and 
the  minimum  energy  requirement  of  a  thermal  or  chemical  pro¬ 
cess  [47,48].  The  hot  composite  curve  represents  all  heat  delivering 
streams  condensed  into  a  single  curve.  Likewise,  the  cold  com¬ 
posite  curve  comprises  all  heat  absorbing  streams.  Combined  in 
the  same  diagram,  they  illustrate  the  overall  heat  transfer  as  it 
would  occur  if  all  heat  was  exchanged  in  a  single  heat  exchanger. 
In  the  present  cases,  all  external  heat  exchangers  of  the  APU  sys¬ 
tem  and  also  the  internal  heat  exchangers  depicted  in  Figs.  1  and  2 
are  taken  into  account.  Fleat  transfer  associated  with  mixing  of  hot 
and  cold  streams,  which  theoretically  could  be  considered  by  intro¬ 
ducing  virtual  heat  exchangers  [49],  is  neglected.  Fig.  6  shows  that 
the  composite  curves  of  alternative  4  extend  significantly  less  into 
the  x-direction  than  those  of  the  base  case.  As  a  consequence  of 


Enthalpy  /  kW 

Fig.  6.  Pinch  diagram  showing  composite  curves  of  base  case  (a)  and  system  alter¬ 
native  4  (b).  Solid  lines:  hot  composite  curves  (heat  delivering  streams);  dashed 
lines:  cold  composite  curves  (heat  absorbing  streams). 
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Table  3 

Power  consumption  (W)  of  blowers/compressors  needed  to  convey  the  indicated 
gas  streams. 


Base  case 

Alt.  1 

Alt.  2 

Alt.  3 

Alt.  4 

Air  supply  to  ATR 

40 

50 

40 

40 

40 

Cathode  air  supply 

418 

481 

519 

24 

29 

Anode  recycle 

- 

- 

36 

- 

36 

Cathode  recycle 

- 

- 

- 

271 

333 

Combustion  air  (burner) 

20 

- 

6 

20 

7 

Cooling  air  (condenser) 

54 

50 

42 

54 

42 

this  decreased  heat  exchange,  the  minimal  required  overall  heat 
exchanging  area  can  be  reduced,  which  is  in  agreement  with  the 
last  row  of  Table  2.  Yet  another  conclusion  can  be  drawn  from  Fig.  6: 
Since  both  composite  curves  are  aligned  to  the  right  without  gen¬ 
erating  pinch  points  or  intersections,  the  thermal  power  required 
to  heat  up  all  cold  streams  can  completely  be  delivered  internally 
from  the  hot  streams. 

If  system  alternative  4  could  be  realized  without  any  para¬ 
sitic  power  losses,  theoretically  an  AC  efficiency  of  35.4%  could  be 
obtained,  which  implies  another  potential  increase  of  1 0.6%.  Thus,  a 
first  step  towards  further  system  optimization  is  the  identification 
of  the  main  contributors  to  system-intrinsic  power  consumption. 
These  are  the  gas  blowers  and  compressors,  the  AC  power  needs 
of  which  are  summarized  in  Table  3.  The  cathode  air  blower  is 
clearly  identified  as  the  main  consumer  in  systems  without  cath¬ 
ode  recycling.  In  contrast,  if  cathode  recirculation  loops  are  present, 
the  blower  inside  the  recycle  loop  contributes  most  to  the  overall 
consumption. 

4.2.  Sensitivity  analysis 


Table  5 

24  factorial  design  -  normalized  factor  levels  of  main  factors  and  binary  interaction 
factors. 


Runj  Factor 


FI 

F2 

FI  2 

F3 

FI  3 

F23 

F4 

FI  4 

F24 

F34 

1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

+1 

+1 

2 

+1 

-1 

-1 

-1 

-1 

+1 

-1 

-1 

+1 

+1 

3 

-1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

4 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

-1 

-1 

+1 

5 

-1 

-1 

+1 

+1 

-1 

-1 

-1 

+1 

+1 

-1 

6 

+1 

-1 

-1 

+1 

+1 

-1 

-1 

-1 

+1 

-1 

7 

-1 

+1 

-1 

+1 

-1 

+1 

-1 

+1 

-1 

-1 

8 

+1 

+1 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

9 

-1 

-1 

+1 

-1 

+1 

+1 

+1 

-1 

-1 

-1 

10 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

-1 

-1 

11 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

-1 

+1 

-1 

12 

+1 

+1 

+1 

-1 

-1 

-1 

+1 

+1 

+1 

-1 

13 

-1 

-1 

+1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

14 

+1 

-1 

-1 

+1 

+1 

-1 

+1 

+1 

-1 

+1 

15 

-1 

+1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

+1 

16 

+1 

+1 

+1 

+1 

+1 

+1 

+1 

+1 

+1 

+1 

-  Y3:  standard  flow  rate  of  air  to  the  cathode; 

-  Y4:  lambda  of  the  afterburner  as  defined  in  Eq.  (18). 

The  factorial  design  (Table  5)  constitutes  a  matrix  with  elements 
Cy  varying  between  -1  and  +1.  According  to  [50],  this  allows  for 
the  calculation  of  quantitative  effects  of  all  factors  i  on  a  response 
variable  Yk : 

N 

£m=qTn  Dcoy*j)  (22) 

j= 1 


In  the  previous  section,  the  system  efficiency  of  an  SOFC-APU 
was  optimized  by  varying  the  process  layout.  Constant  process 
parameters  were  assumed,  with  their  respective  values  adjusted 
to  the  state  of  the  art.  However,  future  technological  development 
will  potentially  yield  improvements  in  most  system  components. 
In  order  to  examine  the  resulting  potential  for  a  further  enhance¬ 
ment  of  system  efficiency,  a  sensitivity  analysis  was  conducted.  As 
indicated  in  Table  4,  four  critical  process  parameters  were  selected. 
Each  of  these  parameters  was  varied  between  two  stages:  the  actual 
value  and  an  assumed  improved  value.  Since  the  aim  was  to  iden¬ 
tify  an  optimized  APU  system,  the  sensitivity  analysis  was  based  on 
process  alternative  4. 

The  schematic  of  the  applied  24  factorial  design  is  outlined  in 
Table  5.  It  contains  four  main  factors  (FI,  F2,  F3  and  F4)  which 
correspond  to  the  process  parameters  introduced  in  Table  4.  The 
larger  and  smaller  values  of  the  two  stage  variations  are  indicated 
by  “+1”  and  “-1”,  respectively.  Furthermore,  all  two-factor  interac¬ 
tions  (FI 2,  FI 3,  etc.)  can  be  assessed  in  a  straightforward  manner, 
which  is  one  major  advantage  of  using  a  factorial  design.From  each 
of  the  16  simulation  results,  the  following  four  response  variables 
were  extracted  for  evaluation: 

-  Y\ :  net  electric  efficiency  of  the  system; 

-  Y2:  lambda  of  the  stack  as  defined  in  Eq.  (10); 


Table  4 

Minimum  and  maximum  values  of  varied  parameters. 


Parameter 

Number 

(Ff) 

Actual  value 
(level  1) 

Improved  value 
(level  2) 

Blower  efficiency 

1 

0.60 

0.75 

Oxygen-to-carbon  ratio  in  ATR 

2 

0.47 

0.45 

Reference  voltage  (mV) 

3 

700 

750 

Stack  fuel  utilization 

4 

0.75 

0.85 

The  sum  in  Eq.  (22),  called  a  contrast,  comprises  all  N=  16  calcu¬ 
lated  values  Ykj.  In  order  to  allow  for  comparison  of  the  effects  on 
different  response  variables,  relative  effects  are  defined: 


^k,i  — 


Lk,i 


CVN)£,V*j 


Ek,i 


(23) 


In  Fig.  7,  the  relative  effects  are  plotted  against  the  four  main 
and  six  interaction  factors. 

First  of  all,  Fig.  7  clearly  shows  that  the  effects  of  the  interaction 
factors  are  small  compared  to  those  of  the  main  factors.  The  effects 
of  the  main  factors  are  discussed  in  detail. 


-  Blower  efficiency  (FI):  Increasing  the  adiabatic  efficiency  of  the 
blowers  from  0.6  to  0.75  improves  the  net  electric  efficiency 
of  the  system  by  approximately  2%.  The  decrease  in  the  power 
dissipation  of  the  adiabatic  blowers  leads  to  lower  blower  out- 
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Fig.  7.  Relative  effects  of  the  main  factors  (cf.  Table  4)  and  interaction  factors  on  net 
electric  efficiency  (Y1 ),  air-to-fuel  ratio  of  the  stack  (Y2),  standard  flow  rate  of  air  to 
the  cathode  (Y3)  and  air-to-fuel  ratio  of  the  burner  (Y4). 
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let  temperatures.  As  a  consequence,  the  amount  of  excess  air 
required  to  cool  the  stack  decreases. 

-  Oxygen-to-carbon  ratio  in  the  ATR  (F2):  More  liquid  fuel  is  con¬ 
verted  by  endothermic  steam  reforming  and  less  by  exothermic 
partial  oxidation  when  n(02)/n(C)  decreases  from  0.47  to  0.45.  As 
a  result,  more  hydrogen  is  produced  per  mole  of  fuel  and  con¬ 
sequently  the  net  electric  efficiency  of  the  whole  system  (Yi ) 
increases.  Furthermore,  reformate  leaves  the  ATR  at  a  temper¬ 
ature,  which  is  lower  by  25  °C,  resulting  in  an  increased  methane 
content.  A  lower  reformate  temperature  and  increased  rates  of 
endothermic  methane  steam  reforming  in  the  SOFC  mean  that 
less  excess  air  has  to  be  supplied  to  the  stack  for  cooling  (cf.  [51  ]). 
In  this  regard,  the  normalized  air-to-fuel  ratio  (Y2)  appears  to  be 
more  affected  than  the  flow  rate  of  air  (Y3).  However,  this  merely 
refers  to  the  fact  that  the  absolute  flow  rate  of  utilizable  fuel  gas 
increases  as  a  result  of  the  higher  rate  of  steam  reforming.  Finally, 
decreasing  the  oxygen-to-carbon  ratio  has  an  unfavorable  effect 
on  the  amount  of  fresh  air  required  for  uncritical  operation  of  the 
afterburner  ( Y4 ).  This  is  caused  by  the  lower  overall  mass  flow  rate 
of  anode  off-gas,  which  simultaneously  contains  a  larger  amount 
of  combustible  compounds. 

-  Reference  voltage  (F3):  Qualitatively,  F3  causes  the  same  effects 
as  FI.  However,  the  absolute  values  of  the  effects  are  almost  five 
times  larger,  which  can  be  ascribed  to  the  large  impact  of  the 
stack  voltage  on  the  response  variables  Y\  to  Y3.  These  are  more 
directly  affected  by  F3  compared  to  FI . 

-  Stack  fuel  utilization  (F4):  Increasing  stack  utilization  from  0.75  to 
0.85  improves  the  net  electric  system  efficiency  by  6%.  However, 
more  air  is  needed  to  maintain  the  specified  operation  tempera¬ 
ture  of  750  °C.  Interestingly,  the  relative  effect  on  the  absolute  air 
flow  is  4.5  times  larger  than  the  effect  on  lambda.  This  is  due  to 
the  definition  of  lambda  in  Eq.  (10),  which  accounts  for  the  uti¬ 
lized  fuel  instead  of  the  total  fuel  flow  rate.  Moreover,  increasing 
the  fuel  utilization  significantly  lowers  the  concentration  of  com¬ 
bustible  compounds  in  the  anode  off-gas.  As  a  consequence,  less 
fresh  air  is  required  to  limit  the  adiabatic  combustion  tempera¬ 
ture  in  the  burner,  resulting  in  a  very  large  effect  of  factor  F4  on 
the  lambda  of  the  afterburner  (Y4). 


standard  flow  rate  of  effective  off-gas  is  as  low  as  25.3  m3  h-1  at 
51 7  °C.  Only  0.47  m2  of  heat  exchanging  area  is  needed.  Thus,  except 
for  the  condensation  temperature,  which  increases  to  a  value  of 
66  °C,  all  process  targets  could  be  improved  by  parameter  optimiza¬ 
tion. 

5.  Conclusion 

A  thermodynamic  model  for  numerical  simulations  of  solid 
oxide  fuel  cell  systems  including  autothermal  reforming  was  pre¬ 
sented.  The  approach  is  based  completely  on  highly  developed 
process  simulation  software,  which  allows  building  and  varying 
complex  system  layouts  with  only  limited  effort. 

A  conceptual  study  was  conducted,  in  which  five  different  pro¬ 
cess  designs  of  an  SOFC-based  auxiliary  power  unit  were  analyzed 
and  compared  to  each  other.  Gains  in  system  efficiency  could  be 
achieved,  however  at  the  cost  of  more  complex  system  designs. 
Compared  to  the  reference  case,  the  high-efficiency  system  con¬ 
tains  a  cathode  and  an  anode  recycle  loop.  In  order  to  evaluate 
further  optimization  potentials,  a  parametric  study  was  conducted, 
including  four  critical  process  parameters  associated  with  different 
system  components.  Using  techniques  from  the  design  and  analysis 
of  experiments,  factor  effects  were  quantified  and  linear  regression 
polynomials  could  be  extracted.  By  combining  the  optimized  pro¬ 
cess  design  and  optimized  process  parameters,  a  net  electric  system 
efficiency  of  38.5%  was  achieved,  which  corresponds  to  an  increase 
of  almost  1 0%  points  compared  to  the  reference  system.  These  theo¬ 
retical  findings  may  serve  as  guidelines  for  the  future  development 
of  prototype  or  even  close-to-production  systems. 

One  of  the  main  obstacles  to  a  further  improvement  of  the  fuel 
cell  system  performance  was  found  in  intrinsic  electric  power  con¬ 
sumption.  The  blowers,  which  deliver  air  to  the  cathode,  were 
identified  as  the  main  consumers.  Thus,  since  air  is  the  primary 
cooling  medium  for  the  fuel  cell  stack,  alternative  stack-cooling 
concepts  are  desirable.  Possible  designs  for  such  alternative  cooling 
concepts  will  be  examined  in  a  forthcoming  study. 
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Only  one  significant  interaction  effect  was  observed:  the  effect 
of  F24  on  the  normalized  air-to-fuel  ratio  of  the  afterburner  (Y4). 
F24  describes  the  interaction  between  the  ATR’s  oxygen-to-carbon 
ratio  (F2)  and  stack  fuel  utilization  (F4);  it  indicates  that  the  effect 
of  F2  on  Y4  is  different  at  different  factor  levels  of  F4.  In  the  present 
case,  however,  this  simply  originates  from  a  process  constraint:  at 
high  fuel  utilization,  i.e.  at  the  upper  level  of  factor  F4,  the  response 
variable  Y4  cannot  adjust  freely  to  F2  but  remains  at  the  minimal 
allowable  value  of  1.1. 

As  stated  at  the  beginning  of  this  section,  the  main  goal  of  the 
sensitivity  analysis  was  to  evaluate  the  potential  for  improving 
system  performance  by  parameter  variation.  For  the  purpose  of 
quantification,  a  first-order  regression  polynomial  was  derived  on 
the  basis  of  the  relative  effects.  Within  the  limits  of  all  /  =  10  factors 
(-1  <  Cj  <  1),  the  response  variables  can  be  calculated  as  follows: 


The  final,  parameter-optimized  system  is  based  on  the  pro¬ 
cess  layout  of  alternative  4  (Fig.  5)  and  is  thus  denoted  alternative 
4*.  Compared  to  alternative  4,  system  fuel  utilization  increases  to 
91.9%  and  intrinsic  electric  power  consumption  decreases  to  only 
0.38  kW.  As  a  consequence,  a  net  electric  efficiency  of  38.5%  is 
reached  for  alternative  4*.  It  exhibits  an  SOFC  lambda  of  2.0,  a  cor¬ 
responding  air  standard  flow  rate  of  19  m3  h-1  and  a  standard  flow 
rate  of  air  to  the  afterburner  of  only  0.9  m3  h-1.  Furthermore,  the 
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